function y = specA2d(N,u,Vrep, ksqrrep)
h = 1 / N;
ud2 = ksqrrep .* fft(fft(u,[],1),[],2)*h^2;
ud2 = ifft(ifft(ud2,[],1),[],2)*N*N;
y = real(Vrep.*u + ud2);
